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Software Development for Raster Change 
Detection using Machine Learning and 
Deep Learning Techniques 


1. Introduction: 


Change detection is the process of identifying differences in the state of an object or 
phenomenon by observing it at different times. Multi-temporal RS data, such as 
satellite imagery and aerial imagery, can provide abundant information to identify 
land use and land cover (LULC) differences in a specific area across a period of time. 
This is very crucial in various applications, such as urban planning, environmental 
monitoring, agriculture investigation, disaster assessment, and map revision. 


2. Change Detection Method: 


Present the data from different sensors including optical RS ( Remote 
Sensing ) data, SAR ( Synthetic Aperture Radar ) data which will be used for 
training and evaluating ML models in change detection method. 


2.1 Image Dataset Information: 


* Jpeg-2000 
* Any Face Image 
* Video Surveillance data 


2.2 Using CNN (Convolutional Neural Networks ): 


A convolutional neural network is a feed-forward neural network that is 
generally used to analyse visual images by processing data with grid-like 
topology. It’s also known as a ConvNet. A convolutional neural network is 
used to detect and classify objects in an image. 


The convolution operation forms the basis of any convolutional neural 
network. Let’s understand the convolution operation using two 
matrices, a and b, of 1 dimension. 


a = [5,3,7,5,9,7] , b=[1,2,3] 


In convolution operation, the arrays are multiplied element-wise, and the 
product is summed to create a new array, which represents a*b. The first 
three elements of the matrix a are multiplied with the elements of matrix b. 


The product is summed to get the result. 


Sum the product 


F a*b| Multiply the arrays 
element wise 


[5, 6, 6] 


a = [5, 3, 2, 
b =[], 2, 3] 


a*b=[17,] 


The next three elements from the matrix a are multiplied by the elements in 
matrix b, and the product is summed up. 


) 
fa*b! Multiply the atley= Sum the product 
element wise 


[S, 6, 6] 7 
a=[5.3,2,5, (3, 4, 15] 22 
b=[1,2,3] 
a*b=[17,22] 


The process continues until the convolution operation is complete. 


In CNN every image is represented in the form of an array pixel values. 


Real Image of the digit 8 Represented in the form Digit 8 represented in the form 
of an array of pixels of O’s and 1’s 


Real Image Represented in the form of Image represented in the 
black and white pixels form of a matrix of numbers 


This is another example to depict how CNN recognises an image. As you 
can see from the above diagram, only those values are lit that have a value 
of 1. 


2.3 Different Layers in Convolutional Neural Network: 


1. Convolution Layer: This is the first step in the process of extracting 
valuable features from an image. A convolution layer has several filters that 
perform the convolution operation. Every image is considered as a matrix of 
pixel values. 


y ~ Image 

/ pixels 

| ~ Filter 
1} 1/0} 0 / 
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Consider the following 5x5 image whose pixel values are either O or 1. 
There’s also a filter matrix with a dimension of 3x3. After Sliding the filter 
matrix over the image and compute the dot product to get the convolved 
feature matrix. 


The stacking of convolutional layer allows a hierarchal decomposition of the 
input. Consider that the filter that operates directly on the row pixel values 
will learn to extract low level features such as lines. By adding more hidden 
layers / more neurons per layer we add more parameters to the model. 
Hence we allow the model to fit more complex function. 


We can Use 2 Api’s Conv2D and Convolution2D to perform convolution. 


2. ReLU Layer ( Rectified Linear Unit ): Next step is to move feature maps 
that are extracted to a ReLU Layer. ReLU performs an element-wise 
operation and sets all the negative pixels to O. It introduces non-linearity to 
the network, and the generated output is a rectified feature map. Below is 
the graph of a ReLU function: 


RZ) = max(O, z) 


The original image is scanned with multiple convolutions and ReLU layers for 
locating the features. 


We use ReLU Api and function from Keras or Tensorflow to perform ReLU 
task. 


3. Pooling Layer: Pooling is a down-sampling operation that reduces the 
dimensionality of the feature map. The rectified feature map now goes 
through a pooling layer to generate a pooled feature map. 


Rectified feature map 


Pooled feature map 


max pooling with 2x2 filters 


and stride 2 


Max(3, 4, 1, 2) = 4 


The pooling layer uses various filters to identify different parts of the image 
like edges, corners, body, feathers, eyes, and beak. 


Here’s how the structure of the convolution neural network looks so far: 
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Max Pooling:- Sharper than original Image 

Average Pooling:- Smoothes the image while keeping the essence of the 
feature in an image 

Minimum Pooling:- Used when image has light background , since min 
pooling will select darker pixels 


We use MaxPooling2D Api and multiprocessing.pool(), max_pool() and 
other similar function to implement pooling. 


These 3 layers, 


1. Convolutional Layer 
2. ReLU Layer 
3. Pooling Layer 


Comes under Feature Learning part of CNN. Now we will talk about our next 
step and next part Classification; 


1. Flattening: \|t takes pooled feature map that is generated in pooling step 
and transform it into a one dimensional vector. Flattening is used to 
convert all the resultant 2-Dimensional arrays from pooled feature maps 
into a single long continuous linear vector. 


Pooled feature map 


Reason to transform pooled feature map into a 1D vector is because this 
vector will now be fed int ANN. This vector will be input of ANN network that 
will be chained onto CNN we will be building. 


We use Flatten Api and convo_2_flat = tf.reshape(), 
tf.keras.layers.Flatten() function to introduce flattening step. 


2. Full Connection: 


1/1}1/0}0 
0/1}/1/1\/0 ; 
1 --—7— Convolution Poolin Flattenin 
O/ 0/1/17 | pe — — 
0;0;1;1)0 
1/1 
2 ae Pooling Layer 
Input Image 


Hidden Layer of the ANN is replaced by a specific type 
of hidden layer called a fully connected layer. A fully connected layers 
neurons are connected to all the neurons in next layer. 


Convolution Layer 


+ 


Fully Connected 
Layer 


Purpose is to detect certain features in image. Each neuron in fully 
connected layer corresponds to specific feature that might be present in an 


image. 


We use full_layer-one = tf.nn.relu(normal_full_layer()) and other functions 
to implement fully connected layer. 


3. Testing Model 


3.1 Using MNIST Dataset: We start basic testing of model by using MNIST 


dataset, which is the set of Gray Scale Handwritten integer and digit images. 


from 
from 
from 
from 
from 
from 
from 
from 
from 
from 


numpy import unique, argmax 


tensorf 
tensorf 
tensorf 
tensorf 
tensorf 
tensorf 
tensorf 
tensorf 
matp lot 


lib 


low. 
low. 
low. 
low. 
low. 
low. 
low. 
low. 


keras 
keras 
keras 


.-datasets.mnist import load_data 


import Sequential 


. Layers import Conv2D 
keras. 
keras. 
keras. 
keras. 
keras. 


layers import MaxPool2D 
layers import Dense 
layers import Flatten 
layers import Dropout 
utils import plot_model 


import pyplot 
import matplotlib.pyplot as plt 
import numpy as np 
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The MNIST handwritten digit classification problem is a standard dataset 
used in computer vision and deep learning. These where some of the libraries 
imported to make the model. Unique & Argmax where imported by numpy. 
Unique Returns the sorted unique elements of an array and Argmax function 
returns indices of the max element of the array in a particular axis. 
Load_data will be further defined in program to load mnist dataset. Next we 
import sequential model using tenserflow.keras, Sequential model is 
appropriate for a plain stack of layers where each layer has exactly one input 
tensor and one output tensor. To use 2D convolutional layer we import 
Conv2D, this layer creates a convolution kernel that is wind with layers input 
which helps produce a tensor of outputs. MaxPool2D library does Max 
Pooling operation for 2D spatial data. Dense imports densely connected 
neural network or fully connected neural network, this layer contains all the 
neurons that are deeply connected within themselves, this means that every 
neuron in the dense layer takes the input from all the other neurons of the 
previous layer. Flatten imports flatten layer which is used to flatten the input 
without affecting the batch size, Flatten layer flattens each batch in the inputs 
to 1-dimension. The Dropout layer randomly sets input units to 0 with a 
frequency of rate at each step during training time, which helps prevent 
overfitting. Inputs not set to O are scaled up by 1/(1 - rate) such that the sum 
over all inputs is unchanged. Pyplot is a collection of command style 
functions that make matplotlib work like MATLAB. Each pyplot function makes 
some change to a figure: e.g., creates a figure, creates a plotting area ina 
figure, plots some lines in a plotting area, decorates the plot with labels, etc. 


# loading the MNIST dataset 

(x_train, y_train), (x_test, y_test) = load_data () 

# reshaping the training and testing data 

x_train = x_train.reshape ((x_train.shape[@], x_train.shape [1], x_train.shape 
x_test = x_test.reshape ((x_test.shape[@], x_test.shape[1], x_test.shape[2], 1) 


x_train = x_train.astype('float32')/255.0 # 255 = pixel values are integer val 
x_test = x_test.astype('float32')/255.@ # changes in 32 bit memory 


In next step we load the MNIST dataset and define the training and testing 
functions. 

x_train: uint8 NumPy array of grayscale image data with shapes (60000, 28, 
28), containing the training data. Pixel values range from 0 to 255. 

y_train: uint8 NumPy array of digit labels (integers in range 0-9) with shape 
(60000,) for the training data. 

x_test: uint8 NumPy array of grayscale image data with shapes (10000, 28, 
28), containing the test data. Pixel values range from 0 to 255. 


AZ 


y_test: uint8 NumPy array of digit labels (integers in range 0-9) with shape 
(10000,) for the test data. 


reshape() function changes the shape of an array, the shape of an array is the 
number of elements in each dimension. By reshaping we can add or remove 
dimensions or change number of elements in each dimension. It defines the 
matrix into 28*28 pixels, and 1 represents the colour channel which is 1D 
grayscale. 


In next step it converts the pixels in range [0,255] to range [0,1],which will 
make images contributes more evenly to the total loss. This involves first 
converting the data type from unsigned integers to floats, then dividing the 
pixel values by the maximum value. Float32 changes it into 32 bit memory. 


fig = plt.figure(figsize= (5,3)) 

for i in range (15): 
ax = fig.add_subplot(2,10, i+1, xticks=[], yticks=[]) 
ax. imshow(np.squeeze(x_train[i]), cmap='gray') 
ax.set_title (y_train [i]) 


5 0 4 1 9 2 1 3 1 4 
HORMMEMEM 


2 o 8 & § 
3} ]3/e} 1 


In this step as we can see, we tried to print the images from first 15 index in 
MNIST dataset. 

Now we will determine the shape of the input images and start adding 
different convolutional neural network to the model which includes Conv2D, 
MaxPool2D, Dense, Flatten and Dropout. 


First we will define model as sequential, Sequential model is appropriate for 
a plain stack of layers where each layer has exactly one input tensor and one 
output tensor. Then we will add a single layer of convolution with a small 
filter size (8,3) and a modest number of filters (82) followed by a max pooling 
layer. The filter maps can then be flattened to provide features to the 
classifier. we Know that we will require an output layer with 10 nodes in order 
to predict the probability distribution of an image belonging to each of the 10 
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classes. This will also require the use of a softmax activation function. 
Between the feature extractor and the output layer, we can add a dense layer 
to interpret the features, in this case with 500 nodes. 


# determine the shape of the input images 
img_shape = x_train.shape [1:] 
print (img_shape) 


(28, 28, 1) 


model = Sequential() 

model.add(Conv2D(32, (3,3), activation='relu', input_shape=img_shape) ) 
model.add(MaxPool2D((2, 2))) 

model.add(Conv2D (48, (3,3), activation='relu') ) 

model.add(MaxPool2D ((2, 2))) 

model. add (Dropout (@.5) ) 

model.add(Flatten() ) 

model.add(Dense(50@, activation='relu')) 

model.add(Dense(10, activation='softmax') ) 

model. summary () 


Model: "sequential_5" 


‘Layer (type) ~—Output Shape =——O—C Parc HC 
“conv2d_4 (Conv2D) —~—~=« (None; 26, 26, 32) «320 
max_pooling2d_2 (MaxPooling (None, 13, 13, 32) 1) 
2D) 
conv2d_5 (Conv2D) (None, 11, 11, 48) 13872 
max_pooling2d_3 (MaxPooling (None, 5, 5, 48) 1) 
2D) 
dropout_1 (Dropout) (None, 5, 5, 48) Q 
flatten_1 (Flatten) (None, 1200) i) 
dense_2 (Dense) (None, 500) 600500 
dense_3 (Dense) (None, 10) 5010 


Total params: 619,702 
Trainable params: 619,702 
Non-trainable params: Q@ 


All layers will use the ReLU activation function, which changes all the 
negative pixels to 0. We will totally use 2 convolutional layer, in second 
convolutional layer we will use 48 modest number of filters for better training. 
After applying drop probability of 0.5, 500 neutrons will be randomly dropped 
in every iteration and flatten layer converts 2D array from pooled feature 
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maps into a single long continuous linear vector. At last using 
model.summary() function to get the output of models. 


model.compile (optimizer='adam', loss='sparse_categorical_crossentropy', metrics 
xX = model. fit(x_train, y_train, epochs=10, batch_size=128, verbose=2, validation 


Epoch 1/10 

422/422 - 50s - loss: @.2365 - accuracy: @.9278 -— val_loss: 0.0561 - val_ac 
Epoch 2/10 

422/422 - 43s - loss: 0.0813 - accuracy: @.9744 - val_loss: @.0398 - val_ac 
Epoch 3/10 

422/422 - 43s - loss: @.@599 - accuracy: @.9809 —- val_loss: 0.0359 - val_ac 
Epoch 4/10 

422/422 - 44s - loss: @.0@472 - accuracy: @.9844 - val_loss: 0.0305 - val_ac 
Epoch 5/10 

422/422 - 43s - loss: @.0409 - accuracy: @.9872 - val_loss: 0.0310 - val_ac 
Epoch 6/10 

422/422 - 44s - loss: @.0356 -— accuracy: @.9881 —- val_loss: 0.0309 - val_ac 
Epoch 7/10 

422/422 - 43s - loss: 0.0324 - accuracy: @.9895 - val_loss: @.0308 - val_ac 
Epoch 8/10 

422/422 - 44s - loss: @.0@278 -— accuracy: @.9909 - val_loss: 0.0261 - val_ac 
Epoch 9/10 

422/422 - 44s - loss: @.@261 - accuracy: @.9914 - val_loss: 90.0256 - val_ac 
Epoch 10/10 

422/422 - 43s - loss: @.@239 - accuracy: @.992@ - val_loss: 90.0211 - val_ac 


loss, accuracy = model.evaluate(x_test, y_test, verbose=0) 
print (f'Accuracy: {accuracy*100}' ) 


Accuracy: 99.27999973297119 


Now its time to compile the model, we used optimizer adaptive movement 
estimation(ADAM) is extension of stochastic gradient descent to update 
network weightage during training, metrics ‘accuracy’ improves accuracy 
and speed of the model. sparse_categorical_crossentropy is used as a loss 
function for multi-class classification model where the output label is 
assigned integer value (0, 1, 2, 3...). 

Next we fit the model using epochs = 10, 128 training instances in each 
batch using verbose as 2 and validation_split = 0.1, for 0.1, 10% data will be 
used in the validation set and 90% data will be used in the test set. And after 
evaluating the model we can se our accuracy of model as 99.279. 


At last we see that we give the index as 5 in program which has value of 
handwritten image 2, and then used model.predict to predict the image. And 
we can see that our model has predicted the image correctly as 2 which 
shows that our first trail model is working perfectly. 
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image = x_train[5] 

# lets display the image which we want to predict 
plt.imshow(np.squeeze(image), cmap='gray') 
plt.show() 


image= image.reshape (1, image.shape [@], image.shape [1], image.shape 


p = model.predict ( [image] ) 
print ('Predicted: {}'.format(argmax (p))) 
1/1 [Ssssssssssssssssssssss========] - Os 109ms/step 


Predicted: 2 


3.2 Using RGB Dataset: Now we will test model by using RGB dataset, 
here we use satellite images of cyclone Biparjoy came in India started from 
east central Arabian Sea Gujarat state region from 04 June 2023 to 20 June 
2023 and travelled through different states like Rajasthan. Here in this 
dataset we will have 2 classes cyclone and monsoon, at the end our model 
will recognise which image had cyclone and which had monsoon. 


So lets start coding to process this RGB dataset, first we will import some 
necessary packages to load preprocessed libraries which will help us. 

First we will import Input, Lambda, Dense and Flatten from 
tensorflow.keras.layers. To take the input from user we use input() function in 
python, it is also used to instantiate a keras tensor. Importing Lambda 
imports the lambda layer which exists so that arbitrary expressions can be 
used as a Layer when constructing Sequential and Functional API models. 
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[2]) 


from tensorflow.keras. layers import Input, Lambda, Dense, Flatten 

from tensorflow.keras.models import Model 

from tensorflow.keras.applications.inception_v3 import InceptionV3 
#from keras.applications.vgg16 import VGG16 

from tensorf low. keras.applications.inception_v3 import preprocess_input 


from tensorflow.keras.preprocessing import image 

from tensorflow.keras.preprocessing.image import ImageDataGenerator, load_img 
from tensorflow.keras.models import Sequential 

import numpy as np 

from glob import glob 


Lambda layers are best suited for simple operations or quick 
experimentation. Then we import Model from tensorflow.keras.models, it 
groups the layers into an object with training/inference features, With the 
"Functional API", where we start from Input, we chain layer calls to specify 
the model's forward pass, and finally we create our model from inputs and 
outputs. A new Functional API model can also be created by using the 
intermediate tensors. This enables you to quickly extract sub-components of 
the model. 


After we import InceptionV3 from tensorflow. keras. applications. 
inception_v3, basically inceptionV3 is an architecture of CNN which is used 
as an image recognition model that has been shown to attain greater than 
78.1% accuracy on the ImageNet dataset. It previously contains 93 layers of 
conv2D, 8 avg pooling and 3 max pooling layers. It instantiates the 
inceptionV3 architecture, this function returns a Keras image classification 
model, optionally loaded with weights pre-trained on ImageNet. Importing 
preprocess_input preprocesses a tensor or numpy array encoding batch of 
images. After importing image iterator is capable of reading images from a 
directory or drive from disk. At last we import glob, In Python, the glob 
module is used to retrieve files/pathnames matching a specified pattern. 


#resize image 

IMAGE_SIZE = [224,224] 

train_path = '/content/drive/MyDrive/Cyclone/training' 
valid_path = '/content/drive/MyDrive/Cyclone/testing' 


# import the inception v3 and add preprocessing layer 

# imagenet weights are used 

inception = InceptionV3(input_shape=IMAGE_SIZE + [3], weights='imagenet', includ 
# false because we dont want to use first and the last layer 


Downloading data from 
87910968/87910968 [ 
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Now we will resize the image in dataset to [224, 224], because it reduces the 
size of large files which makes it easier to process. Provide path of training 
and testing dataset for which we have to resize and define it as train_path 
and valid_path. Then we will import inceptionV3 model and add 
preprocessing layer where imagenet weights are used. While loading a given 
model, the “include_top” argument can be set to False, in which case the 
fully-connected output layers of the model used to make predictions is not 
loaded, allowing a new output layer to be added and trained. Optional shape 
tuple input_shape can only be specified when include_top is false, it should 
have exactly 3 input channels (224, 224, 3) where 224*224 specifies width 
and height which is in a matrix of 224*224 and 3 defines the colour channel 
or we can say 3D array, width and height should be not smaller than 75. 


#for not training the existing weights 
for layer in inception. layers: 
layer.trainable = False 


#useful for getting number of output classes 
folders = glob('/content/drive/MyDrive/Cyclone/training/x' ) 
folders 


['/content/drive/MyDrive/Cyclone/training/cyclone', 
'/content/drive/MyDrive/Cyclone/training/monsoon' ] 


Putting layer.trainable to false moves all the layers weight from trainable to 
non trainable which is freezing the layer, the state of the frozen layer won’t be 
updated during training. Then we will use glob to Know number of existing 
classes in training dataset. 


x = Flatten()((inception. output) ) 

prediction = Dense(len(folders), activation='softmax') (x) # softmax for multig 
# creating a model object 

model = Model(inputs=inception. input, outputs=prediction) 

model. summary ( ) 


batch_normalization_87 (Ba (None, 5, 5, 384) ['conv2 
tchNormalization) 


batch_normalization_88 (Ba (None, 5, 5, 384) ['conv2 
tchNormalization) 


batch_normalization_91 (Ba (None, 384) ['conv2 
tchNormalization) 


batch_normalization_92 (Ba (None, 384) ['conv2 
tchNormalization) 


conv2d_93 (Conv2D) (None, 192) 393216 ['avera 


batch_normalization_85 (Ba (None, 320) 960 ['conv2 


Now we will apply flatten to the output of inception layer, flattening is 
converting the data into a 1D array for inputting it to the next layer, we flatten 
the output of the convolutional layers to create a single long feature vector, 
here we can add more layers if we want. After we create dense layer which is 
a fully connected neural network with activation function softmax, it will 
convert linear output into probabilistic. Then we derived the model, and 
“model = Model(inputs=inception.input, outputs=prediction)” this is the 
model which we are going to train and model.summary shows up the every 
layer detail of model in output. 


, , 
# adpti sio r oda weightage during training 
# using or th 

train_datagen = ImageDataGenerator(rescale = 1./255,shear_range = Q.2,zo0om_range = @.2,horizontal_flip = True) 

test_datagen = ImageDataGenerator(rescale = 1./255) 

validation_datagen = ImageDataGenerator (rescale= 1./255) 

# convert the pixels in range [0,255] to range [0,1],will make images contributes more evenly to the total loss 


from google.cola rt drive 
drive.mount( '/cor en 


Drive already mounted at /content/drive; to attempt to forcibly remount, call drive.mount("/content/drive", force_remount=True). 


training_set = train_datagen.flow_from_directory(' 


test_set = test_datagen. flow_from_directory( '/c: 
tar 
validation_set = validation_datagen. flow_from_directory(' 
224,224),batch_size = 18,class_mode = al') # categorical=2d 1-hot encoding label 
r = model. fit_generator( training_set, validation_data = val _set, 
epochs=500, verbose=2, steps_per_epoch=len(training_set) , validation_steps=len(validation_set) ) 


1/1 - 2s - loss: @.0000e+00 - accuracy: 1.0000 - val_loss: 0.@000e+@@ - val_accuracy: 1.0000 - 2s/epoch - 2s/step 
Epoch 472/500 
1/1 - 2s - loss: @.0000e+00 — accuracy: 1.0000 - val_loss: @.@000e+00 - val_accuracy: 1.0000 - 2s/epoch — 2s/step 
Epoch 473/500 
1/1 - 1s - loss: @.0000e+00 — accuracy: 1.0000 - val_loss: 0.@000e+00 - val_accuracy: 1.0000 - 1s/epoch — 1s/step 
Epoch 474/500 
1/1 - 1s - loss: @.000@e+00 - accuracy: 1.0000 - val_loss: 0.000@e+00 - val_accuracy: 1.0000 - 1s/epoch - 1s/step 
Epoch 475/500 
/1 - 1s - loss: @.0000e+0@ — accuracy: 1.0000 - val_loss: 0.0000e+00 - val_accuracy: 1.0000 - 1s/epoch — 1s/step 
Epoch 476/500 
/1 - 2s - loss: @.0@00e+0@ - accuracy: 1.0000 - val_loss: @.000@e+00 - val_accuracy: 1.00@@ - 2s/epoch - 2s/step 
Epoch 477/500 
1/1 - 2s - loss: @.0000e+00 — accuracy: 1.0000 - val_loss: 0.@000e+0@ - val_accuracy: 1.0000 - 2s/epoch — 2s/step 
Epoch 478/500 
1/1 - 2s - loss: @.0000e+00 — accuracy: 1.0000 - val_loss: 0.@000e+00 - val_accuracy: 1.0000 - 2s/epoch — 2s/step 


Now we will compile the model with categorical_crossentropy as a loss 
function, Adaptive Movement Estimation (ADAM) as optimizer which is 
extension of stochastic gradient descent to update network weightage 
during training, and metrics as accuracy, it calculates how often reductions 
equal labels. This metric creates two local variables, total and count that are 
used to compute the frequency and This frequency is ultimately returned as 
binary accuracy. 


Then we will use ImageDataGenerator to import the images from dataset, 
converting the pixels in range [0, 255] to range [0, 1] will make images 
contribute more evenly to the total loss, shear_range = 0.2 means shearing / 
cutting the image by 20%, zoom_range = 0.2 defines zoom in and zoom out 
by 20% and horizontal_flip = true means it will take both rows and columns 
of such a matrix and flips them horizontally along the Y-axis. We will do the 
for different classes and name it as train_datagen, test_datagen, 
validation_datagen. Now because we are using Google colab platform for 
trial we will import drive from google.colab and mount the drive to the server. 
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Now lets create train, test and validation data generator with train, test and 
validation directory as input. flow_from_directory method allows us to read 
the images directly from the directory and augment them while the neural 
network model is learning on the training data. As we know target_size 
defines the size of image we need. Here we need [224, 224] size image so 
that’s the input there, batch_size defines to how pics will be there in every 
batch, and class_mode categorical determines the type of label arrays that 
are returned, categorical is a2D one hot encoded labels. Then at last we fit 
the model, where training_set is the set of data on which actual learning is 
taking place, validation_set helps to improve the model performance by fine- 
tuning the model after each epoch, epochs = 500 defines that the learning 
algorithm will work 500 times through the entire training set. By verbose = 2 
mode verbosity is defined as single line, steps_per_epoch = len(training_set) 
allows us to handle any situation where the number of samples in our epoch 
is different, like here we gave the parameter as training_set and at last 
validation_steps = len(validation_set) tells us about a single backpropagation 
performed on the validation data. we can see our model has an overall 
accuracy of 1.0000 with no loss. 


from matplotlib import pyplot 
import matplotlib.pyplot as plt 
import numpy as np 


img = image. load_img("/content/drive/MyDrive/Cyclone/training/cyclone/1. jpeg") 
plt. imshow( img) 


<matplotlib.image.AxesImage at 0x7fcb58306350> 


import cv2 
cv2.imread("/content/drive/MyDrive/Cyclone/testing") 


training_set.class_indices 


{'cyclone': @, 'monsoon': 1} 
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Now we import pyplot from matplotlib, and then check if image path is 
working perfectly by printing one of the image from library. Then we import 
cv2 which is the module import name for open cv python. Now cv2.imread() 
method loads the images from specified file. 


Then we use training_set.class_indices to get the prediction in binary form, 
As the prediction will be in a binary form, we will be receiving either a O or 1, 
which will represent a cyclone or a monsoon respectively. 


© import os 
from numpy import argmax 
dir_path = ("/content/drive/MyDrive/Cyclone/testing") 
for i in os.listdir(dir_path): 
img = image. load_img(dir_path+'//'+i, target_size = (224,224)) 
plt. imshow( img) 
plt.show() 


X = image. img_to_array(img) 

X = np.expand_dims(X,axis =0) 

images = np.vstack( [X] ) 

p = model.predict(images) 

print ('Predicted: {}'.format(argmax(p) ) ) 


Predicted: 


First we import os, which provides functions for creating and removing a 
directory, fetching its contents, changing and identifying the current 
directory, etc. It directly interacts with the underlying operating system. Then 
from numpy we import argmax, which returns indices of the max element of 
the array in a particular axis. 
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Then we define our directory path for testing, define a for loop and load 
images by defining target size as [224, 224]. Function img_to_array converts 
the PIL image instance to numpy array, with the help of np.expand_dims() 
method we can get the expanded dimensions of an array and axis = 0 show 
that we are talking about multi-dimensional arrays, axis 0 is the axis that runs 
downward down the rows. 

np.vstack() function is used here to stack the sequence of input arrays 
vertically to make a single array. Then model.predict(), predicts the label of a 
new set of data. In print (‘Predicted: {}’.format(argmax(p))) the argmax 
function determine the index position into an array which contains the max 
value and, the highest probability prediction. We can see in above fig that 
our model predicts the cyclone correctly when input image for testing is 
given. Which means model is working absolutely perfect. 


Predicted: 1 


Model perfectly predicts monsoons 
22 


3.3 Predicting Stages of Cyclones 


Now as we got success in processing RGB dataset, lets do a slight change 
in our model and follow our main intution of working. Our main intuition is to 
predict different stages of cyclone just by providing an raster image data to 
the model and then showing change detection in two raster images of same 
place and situation with different time intervals. 


a @& Ts] > [Users/kunalpathako826 
[Users/kunalpathak9826/untitled3.py ; ; GER 


SX change.py  untitied2.py | untitiedsipy 


#for not training the existing weights 
for layer in inception. layers: 
layer.trainable = False 


#useful for getting number of output classes 
folders = glob( '/Users/kunalpathak9626/Desktop/Cyclone/training/#') 
folders 


x = Flatten()((inception-output)) 

Prediction = Dense(len(folders), activation='softmax’) (x) # softnax for multiple categories 
# creating a model object 

model = Model (inputs=inception.input, outputs=prediction) 

model. sunmary() 

#92 layers of conv2D/ 8 avg_pooling/ 3 max_pooling/ 1 fully connected layer 


model. compile(loss= 'categorical_crossentropy’, optimizer=‘adam', metrics={ 'accuracy']) 

# adptive movement estimation(ADAM) is extension of stochastic gradient descent] to update network weightage during tr 

# using image data generator to import the images from dataset 

train_datagen = InageDataGenerator(rescale = 1./255,shear_range = 0.2,zoom_range = 0.2,horizontal_flip = True) 
datagen = ImageDataGenerator(rescale = 1./255) 

validation_datagen = ImageDataGenerator (rescale= 1./255) 

# convert the pixels in range [0,255] to range [@,1],will make images contributes more evenly to the total loss 


training_set = train_datagen. flow_from_directory( '/Users/kunalpathak9826/Desktop/Cyclone/training', 
target_size = (224, 224),batch_size = 18,class_mode = ‘categorical') 
test_set = test_datagen. flow_from_directory( ‘/Users/kunalpathak9626/Desktop/Cydtone/testing * 
target_size = (224, 224) ,batch_size = 18,class_mode = ‘categorical') 
validation_set = validation_datagen.flow_from_directory( '/Users/kunalpathak9826/Desktop/Cyclone/validation 
rget_size = (224,224),batch_size = 18,class_mode = ‘catego: 
r = model. fit_generator( training set, validation_data = validation_set, Help| Variable Explorer | Plots) Files 
epochs=50, verbose=2, steps_per_epoch=len(training_set) ,validation_steps=len(validation_set) 


a 
from matplotlib import pyplot * Console 1/A 


d Speer eres, 
eer ney eecae a js - loss: 0.0168 - accuracy: 1.0000 - val_loss: 0.0068 - val_accuracy: 1.0000 - 8s/epoch - 4s/step 


2/2 
ing = image. load_ing( “/Users/kunalpathak9826/Desktop/Cyctone/testing/stage1(2). ene 


js - loss: 0.0115 ~ accuracy: 1.0000 ~ val_loss: @.@065 - val_accuracy: 1.0000 - 8s/epoch - 4s/step 
plt.imshow (ing) ipocn 45/50 


Keorttcd 2/2 = 8s = loss: 0.0113 - accuracy: 1.0000 - val_loss: 0.0063 - val_accuracy: 1.0000 - &s/epoch ~ 4s/step 
aineget clack aotearoa oe ye ie cee! 2/2 Be = toss: 0.0091 — accuracy: 1.0000 — val_loss: 0.0061 — val_accuracy: 1.0000 ~ 8s/epoch ~ 4s/step 
dain. path = (“/content/drive/MyOrive/Cyclone/testing") a oss: 0093 - accuracy: 1.0000 - val_loss: 0.0059 - val_accuracy: 1.0000 - Bs/epoch - 4s/step 
ein 2 sge toa iat '/Users/kunalpathak9826/Desktop/Cyclone/testing/stage2. jpeg’, target_size = (224,224) 2 = 8s = loss: 0.0088 - accuracy: 1.0000 - val_loss: 0.0057 - val_accuracy: 1.0000 - 8s/epoch - 4s/step 


pte. dnshow ing Epoch 49/50 
plt.show() 2/2 = 8s = loss: 0.0104 - accuracy: 1.0000 - val_loss: @.0055 - val_accuracy: 1000 - 8s/epoch - 4s/step 


Epoch 50/50 

Snoge:4ag_ta_acray lagi 2/2 = 85 = loss: 0.0097 - accuracy: 1.0000 - val_loss: 0.0053 - val_accuracy: 1.000 - 8s/epoch - 4s/step 

np.expand_dins(X, axis=0) 

val = model. predict (X) 

#images = np.vstack( [X]) 

#p = model.predict (images) 

#print ('Predicted: *, format (argmax' 

TF valeny() a ere Figures are displayed in the Plots pane by default. To make them also appear inline in the conso' 
hes lute inline plotting" under the options menu of Plot: 


Xs 
x 


= 05 209ms/step 


This version of python seems to be incorrectly compiled 
(internal generated filenames are not absolute). 

This may make the debugger miss breakpoints. 

Related bug: http://bugs. python. org/issuel666807 


(NO Cyctone") (1 
a = E E IPython Console, History 


conda: base (Python 3.11.3) ‘% Completions: conda(base) ~ LSP:Python Line99,Col1 UTF-8 LF RW Mem75% 


This image shows model working perfectly in an offline setup, using Spyder 
for the development. Model correctly predicts the given image as it is Stage 
2 cyclone as we Can see in image. 


Now lets understand the code and working of the model briefly using 
example of Google Colab Laboratory. 
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r > ) from matplotlib import pyplot 
import matplotlib.pyplot as plt 
import numpy as np 


img = image. load_img("/content/drive/MyDrive/Cyclone/training/Stage 1/paul_tmo_2010088_l1rg.jpg") 
plt. imshow( img) 


<matplotlib.image.AxesImage at 0x7£0d15377fa0> 


1000 2000 3000 4000 5000 


First we import pyplot from matplotlidb and then matplotlib.pyplot as pit, 
pyplot is a collection of command style functions that make matplotlib work 
like MATLAB. Each pyplot function makes some change to a figure: e.g., 
creates a figure, creates a plotting area in a figure, plots some lines in a 
plotting area, decorates the plot with labels, etc. Then we import numpy as 
np, NumPy is a general-purpose array-processing package. It provides a 
high-performance multidimensional array object and tools for working with 
these arrays. It is the fundamental package for scientific computing with 
Python. 


Now we define a function img and try to load a certain jpg image using 
image.load_img() method and then put the path of the image in the bracket 
to load the image. Then we will use plt.imshow(img) method which will 
display the certain image in the console. And we can see that we 
successfully imported the image which means we are in a right direction. 
While working in the google colaboratory we should first focus on mounting 
the google drive successfully so that we can retrieve the data from drive. 


We can do this by using drive.mount() method and passing the certain path 
with the method. 
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import cv2 
cv2. imread("/content/drive/MyDrive/Cyclone/testing") 


training_set.class_indices 


{'Stage 1': 


Qe, 
1, 
A op 
3, 
4, 


‘no_cyclone': 5} 


import os 
from numpy import argmax 


for i in os.listdir(): 
img = image. load_img('/content/drive/MyDrive/Cyclone/testing/stage2.jpeg', target_size = (224,224)) 
plt. imshow( img) 
plt.show() 


image. img_to_array (img) 
np.expand_dims(X,axis=0) 


Xe= 
X= 
val = model.predict(X) 


if val.any() == 0: 
print ("Stage 1") 
elif val.any() == 1: 
print ("Stage 2") 
elif val.any() == 2: 
print ("Stage 3") 
elif val.any() == 3: 
print ("Stage 4") 
elif val.any() == 4: 
print ("Stage 5") 
else: 
print ("N@ Cyclone") 


Next we import cv2, this is an old interface in old OpenCV versions which 
was named as cv. cv2.imread() method loads an image from the specified 
file, but here as we gave the path of a folder which contains 6 different 
folders where each folder contains many RGB images, so it reads that 
folders and returns back without showing any output. 


Training_set.class_indices attribute is used to see the dictionary containing 
the mapping from the class name to class indices. We can see as we pass 
the attribute it gives name of different class stored in training set with their 
indices from 0 to 5. 


Now first we import os, which provides functions for creating and removing a 
directory, fetching its contents, changing and identifying the current 
directory, etc. It directly interacts with the underlying operating system. Then 
from numpy we import argmax, which returns indices of the max element of 
the array in a particular axis. 

Then we define our directory path for testing, define a for loop and load 
images by defining target size as [224, 224]. Function img_to_array converts 
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the PIL image instance to numpy array, with the help of np.expand_dims() 
method we can get the expanded dimensions of an array and axis = 0 show 
that we are talking about multi-dimensional arrays, axis O is the axis that runs 
downward down the rows. 


Then we define a function val = model.predict(X), where model.predict() is a 
function that generates output predictions for the input samples. By default, 
it returns the probabilities of each class for the given input. 


Now we use if and else loop to determine our conditions to how to print our 
predictions as we can see in above image. These conditions determines the 
severity of cyclone, model will determine if cyclone is stage 1, stage 2, stage 
3, stage 4, stage 5 or no cyclone to the given input image accurately. 


if val.any() == @: 
print ("Stage 1") 
elif val.any() == 1: 
print ("Stage 2") 
elif val.any() == 2: 
print ("Stage 3") 
elif val.any() == 3: 
print ("Stage 4") 
elif val.any() == 4: 
print ("Stage 5") 
else: 
print ("NO Cyclone") 


In above image we can clearly see that model predicts the severity of image 
accurately as Stage 2. So our first motive to predict the severity of cyclone is 
successful. 
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4. Raster Change Detection 


Change detection is one of the fundamental applications in imagery and 
remote sensing. It is the comparison of multiple raster datasets, typically 
collected for one area at different times, to determine the type, magnitude, 
and location of change. Change detection is a fundamental requirement of 
remote sensing and raster processing. It is the process of comparing 
multiple rasters or images captured for one location at multiple times to 
identify pixels that have changed due to long-term, seasonal, or abrupt 
changes. ArcGIS Pro allows you to compare rasters and identify the 
magnitude and type of change using tools and raster functions available with 
the Image Analyst extension. 


Here we are using Categorical Raster Data, Categorical raster data is raster 
data where each pixel has a value that is representative of a class or 
category. It is sometimes referred to as discrete data, thematic data, or 
discontinuous data, and it is most often used in GIS to represent land cover, 
land use, or other zonal information such as risk level. The purpose of 
comparing categorical rasters is to identify areas that have changed from 
one class to another class over a period of time. 


27 


4.1 Change Detection Tutorial: 


So, here we begin to write the code to show change detection in two 
different raster images, first task is to install all the required files that we will 
be using here in our program. Here | didn’t had the rasterio package so | had 
to install it using pip install rasterio, you can install packages as per your 
system dependencies. 


pip install rasterio 


Collecting rasterio 
Downloading rasterio-1.3.8-cp310-—cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (21.3 MB) 
213/213 MB eta 0:00:00 


Collecting affine (from rasterio) 

Downloading affine-2.4.@-py3-none-any.whl (15 kB) 
Requirement already satisfied: attrs in /usr/local/lib/python3.10/dist-packages (from rasterio) (23.1.0) 
Requirement already satisfied: certifi in /usr/local/lib/python3.10/dist-packages (from rasterio) (2023.7.22) 
Requirement already satisfied: click>=4.@ in /usr/local/lib/python3.1@/dist-packages (from rasterio) (8.1.6) 
Requirement already satisfied: cligj>=0.5 in /usr/local/lib/python3.10/dist-packages (from rasterio) (0.7.2) 


Requirement already satisfied: numpy>=1.18 in /usr/local/lib/python3.10/dist-packages (from rasterio) (1.22.4) 
Collecting snuggs>=1.4.1 (from rasterio) 

Downloading snuggs-1.4.7-py3-none-any.whl (5.4 kB) 
Requirement already satisfied: click-plugins in /usr/local/lib/python3.10/dist-packages (from rasterio) (1.1.1) 
Requirement already satisfied: setuptools in /usr/local/lib/python3.10/dist-packages (from rasterio) (67.7.2) 
Requirement already satisfied: pyparsing>=2.1.6 in /usr/local/lib/python3.10/dist-packages (from snuggs>=1.4.1— 
Installing collected packages: snuggs, affine, rasterio 
Successfully installed affine-2.4.0 rasterio-1.3.8 snuggs-1.4.7 


Basically Rasterio is a GDAL and Numpy-based Python library designed to 
make us work with geospatial raster data more productive, more fun — more 
Zen. It's a new open source project from the satellite team at Mapbox. 
Rasterio is Python software, not GIS software. 


#Import all the modules we'll need later on 

%Smatplotlib inline 

import numpy as np # Numpy is Python's main numeric processind library 
import datetime # Allows handling of dates 

import matplotlib.pyplot as plt # Add the plotting libraies 


import pyproj 

import rasterio 

import rasterio. features 
import folium 

import ee 

from IPython import display 


# The paths to the Australian Seasonal Landsat Mosaics on the Queensland RDSI Node 
refDataPath = '/vsicurl/http://qld.auscover.org.au/public/data/ landsat/surface_ref lectance/aus/180lre_aus_m201609201 


# Open the surface reflectance data using rasterio (Could also use GDAL or RIOS or rsgislib etc) 
refDataSet = rasterio.open(refDataPath) 


# The dataset object contains metadata about the raster data 
print("Bands:\t ", refDataSet.count) 

print("Height:\t ", refDataSet.height) 

print("Width:\t ", refDataSet.width) 

print("Number of MegaPixels:\t ", refDataSet.heightxrefDataSet.width/1e6,Q) 
print("BoundingBox: ", refDataSet. bounds) 

print("CRS:\t ", refDataSet.crs) 
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So, now we will import all the modules we will need later on. First we import 
Ymatplotlib inline, it enables inline plotting where ploted graphic appear in 
our notebook. Suppose we want to draw a simple line graph using the 
following code. The code works fine, but it doesn't show the line graph inline 
with the code. 


Then we import numpy as np, import datetime, so that our program can 
work with date and times, it is a built in python module. Then import 
matplotlib.pyplot as plt, import pyproj which performs cartographic 
transformations. It converts from longitude, latitude to native map projection x, 
y coordinates and vice versa. After we import rasterio, import 
rasterio.features, import folium which helps us create several types of 
Leaflet maps. By default, Folium creates a map in a separate HTML file. 
Since Folium results are interactive, this library is very useful for dashboard 
building. And now we import ee which allows us to interact with google earth 
engine using the python programming language. At last we import display 
from |IPython, It will provide an IPython terminal and web-based (Notebook) 
platform for Python computing. It has more advanced features than the 
Python standard interpreter and will quickly execute a single line of Python 
code. 


refDataPath = '/v 


refDataSet = rasterio.open(refDataPath 


print("Bands:\t ", refDataSet.count 
print("Height:\t ", refDataSet.height 
print("Width:\ , refDataSet.width 


print ("Number legaPixels ', refDataSet.height*refDataSet.width/1e6, 0 
print("BoundingBox , refDataSet.bounds 
print("CRS: ', refDataSet.crs 


Bands: 6 

Height: 135159 

Width: 141481 

Number of MegaPixels: 19122.430479 0 

BoundingBox: BoundingBox(left=-1944645.@, bottom=-4910195.0, right=2299785.0, top=-855425.0) 
CRS: EPSG:3577 


Now moving to next part, refDataPath function refers to the paths to the 
Australian Seasonal Landsat Mosaics on the Queensland RDSI Node, then 
refDataSet = rasterio.open(refDataPath) opens the — surface reflectance 
data using rasterio (Could also use GDAL or RIOS or rsgislib etc). 


Using print we print metadata such as bands, height, width, number of mega 
pixels, bounding box and CRS about the raster data as we can see in above 
image. 
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# This takes a subset over an AOI with the bounding box in EPSG:3577 (Australian| Albers) 
# Brisbane 

topLeft = (2030000,-3140000) 

bottomRight = (2070000,-3160000) 


# Compute the image coordinates from the real world coordinates using an inverse) transform 
yl, x1 = «refDataSet.transform * topLeft 
y2, x2 = «refDataSet.transform * bottomRight 


# Read in the SWIR, NIR and Red bands within the specified bounding box 
refDataSubset = refDataSet.read([5,4,3], window=((int(x1),int(x2)), (int(y1),int(y2)))) 


# Print some inf 
print( ‘Dat ape: {@}, dty 1}'.format(refDataSubset.shape, refDataSubset.dtype) ) 


# Stretch and scale the data from @ to 255 as an unsigned 8 bit integer for displlay 
refDataSubsetScaled=np.clip(refDataSubset / 4000.0 * 255.0, @, 255).astype('uint8') 


# Plot the image using matplotlib 
plt. figure(figsize=(16,10) ) 


# Rearrange the axes so that it works with matplotlib that likes BIP not BSQ data 
plt.imshow(np. rollaxis(refDataSubsetScaled,@,3) ) 


# Add a title to 


Data extract shap (3, 667, 1333), dtype: int16 
Text(0.5, 1.0, ‘Landsat Surface Reflectance (Bands 5,4,3)') 


Using topLeft = (2030000,-3140000), bottomRight = (2070000,-3160000), 
we takes subset over an AOI with the bounding box in EPSG: 3577 
( Australian Albers ) Brisbane. Then using functions y1, x1 = 
~refDataSet.transform * topLeft, y2, x2 = ~refDataSet.transform * 
bottomRight we compute the image coordinates from the real world 
coordinates using an inverse transform. Function refDataSubset reads the 
SWIR, NIR and Red bands within the specified bounding box, then next print 
functions prints the information of Data Extract Shape, Data Type and 
Landsat Surface Reflectance Bands on the subset. 


Then function refDataSubsetScaled stretch and scale the data from 0 to 255 
as an unsigned 8 bit integer for display. After we plot the image using 
matplotlio and provide the size of fig we need. plt.imshow rearrange the 


Landsat Surface Reflectance (Bands 5,4,3) 
7 ei gi es rf oe PITY 5 / F t 


axes so that it works with matplotlib that likes BIP not BSQ data and at last 
pit.title adds a title to the plot. Above image shows the output of particular 
code, which displays Landsat Surface Reflectance image with Bands 5, 4, 3. 


# Function to convert geographical to projected coordinates 
def 112albers(x,y): 

inProj = pyproj.Proj(init='epsg:4326') 

outProj = pyproj.Proj(init='epsg:3577') 

return pyproj.transform(inProj,outProj,x,y) 


# Utility Query Mapping function 
# Takes a topLeft and bottomRight Tuple 
def mapQuery(topLeft, bottomRight) : 
map_hybrid = folium.Map( 
location=[np.mean( [topLeft[1],bottomRight[1]]), np.mean([topLeft[@],bottomRight [0]])1, 
zoom_start=12, 
tiles="_http://mt1. google. com/vt/lyrs=y&z={z}&x={x}&y={y}", 
attr="Google") 
map_hybrid.add_children(folium. features.PolyLine( 
locations=[ 
(topLeft[1],topLeft[0]), 
(bottomRight [1], topLeft[@]), 
(bottomRight[1] ,bottomRight[0]), 
(topLeft [1], bottomRight[@]), 
(topLeft[1],topLeft[0])], 
color='red', opacity=0.5) ) 
map_hybrid.add_children(folium. features.LatLngPopup() ) 
return map_hybrid 


mapQuery ((150.3727,-26.9819) , (150.4963, -27.0767) ) 


Now we define a function with name Il2albers with attributes x and y. This 
function overall converts given geographical coordinates to projected 
coordinates. Inside ll2albers we define 2 functions inProj & outProj, inProj 
defines EPSG:4326 and outProj defines EPSG:3577 which is a latitude/ 
longitude coordinate system based on the Earth's center of mass, used by 
the Global Positioning System among others. And return fuction returns both 
the function. One of the constants of spatial data is that it's often projected. In 
this case, the projection used is Australian Albers or EPSG:3577. It's 
convenient to specify extract areas in Latitude/Longitude so we define a 
function to do the transform for us. 


In next step we define utility query mapping function as mapQuery which 
takes 2 tuples topLeft and bottomRight. It's useful to define a map to get 
coordinates from, and to plot up bounding boxes to check what we're 
looking at, this map basically used the google hybrid layer. In map_hybrid 
folium.map defines location, zoom start, attribute. Then we add folium 
features in map_ hybrid using map_hybrid.add_children function and at last 
we call the function map_hybrid using return. Then as you can see we 
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passed the coordinates in mapQuery, here we extract the data over an area 
where there is active CSG activity. We can see extracted area in red box in 
below image: 
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# This is a subset west of Tara where Coal Seam Gas development has been happening 
topLeft = 112albers(150.3727,-26.9819) 
bottomRight = 112albers(150.4963,-27.0767) 


# The paths to the 2007 Australian Seasonal Landsat Mosaics on the Queensland RDSI Node 
date1DataPath = '/vsicurl/http://qld.auscover.org.au/public/data/landsat/surface_reflectance/aus/lztmre_aus_m20070324 


# The paths to the 2017 Australian Seasonal Landsat Mosaics on the Queensland RDSI Node 
date2DataPath = '/vsicurl/http://qld.auscover.org.au/public/data/landsat/surface_reflectance/aus/180lre_aus_m2017032q 


# Compute the image coordinates from the real world coordinates using an inverse transform 
y1, x1 ~refDataSet.transform * topLeft 
Wan 3e2 srefDataSet.transform * bottomRight 


# Read in the SWIR, NIR and Red bands within the specified bounding box 
datel1DataSubset = rasterio.open(date1DataPath) .read([5,4,3], window=((int(x1),int(x2)), (int(y1),int(y2)))) 
date2DataSubset = rasterio.open(date2DataPath).read([5,4,3], window=((int(x1),int(x2)), (int(y1),int(y2)))) 


# Setup the plot 
f, (ax1, ax2) = plt.subplots(1, 2, figsize=(16,10) ) 


# Rearrange and scale the data 
ax1.imshow(np.rollaxis(np.clip(datelDataSubset / 5000.@ * 255.0, @, 255).astype('uint8'),@,3), aspect=1) 
ax2. imshow(np.rollaxis(np.clip(date2DataSubset / 5000.@ * 255.0, 0, 255).astype('uint8'),®,3), aspect=1) 


# Add some titles 
ax1,set_title('Autumn 2007 Landsat', fontsize=18) 
ax2.set_title('Autumn 2017 Landsat', fontsize=18) 


topLeft = Il2albers(150.3727, -26.9819) & bottomRight = Il2albers(150.4963, 
-27.0767) are a subset west of Tara where coal seam gas development has 
been happening. date1DataPath function represents the paths to the 2007 
Australian Seasonal Landsat Mosaic on the Queensland RDSI node. 
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Whereas date2DataPath function represents the paths to the 2017 Australian 
Seasonal Landsat Mosaics on the Queensland RDSI node. Our main 
purpose of defining function y1,x1 = ~refDataSet.transform * topLeft , y2,x2 
= ~refDataSet.transform * bottomRight is to compute the image coordinates 
from the real world coordinates using an inverse transform. 


date1DataSubset & date2DataSubset reads in the SWIR, NIR and Red bands 
within the specified bounding box. F, (ax1, ax2) setups the plot, whereas 
ax1.imshow & ax2.imshow rearranges and scales the data, ax1.set_title & 
ax2.set_title adds title to our output image. 


Autumn 2017 Landsat 
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Above image shows the output of 2007 Landsat and 2017 Landsat 
of same place and same location. 


Here we compute NDVI from the extracted subsets, the normalised 
differencing vegetation index is widely used as an indicator of vegetation 
vigour and density. Everyone uses NDVI as “ Hello World ” remote sensing 
example so | will as well, note that these data uses 32767 as a NoData value 
So we need to mask these out from the calculation. 
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# Allow division by zero 
np.seterr(divide='ignore', invalid='ignore') 


# Create a function to compute the ndvi from the red and NIR bands 
def makeNDVI(b4,b3): 
# Calculate NDVI 
ndvi = (b4.astype(float) - b3.astype(float)) / (b4.astype(float) + b3.astype( float) ) 
# Make the NoData values NaN so they get masked from the analysis 
ndvilb4 == 32767] = np.nan 
return ndvi 


# Use the function to produce NDVI images 

# Watch the band ordering — remember we've read in bands 5,4,3 
date1NDVI = makeNDVI(date1DataSubset [1] ,date1DataSubset [2] ) 
date2NDVI = makeNDVI(date2DataSubset [1] ,date2DataSubset [2] ) 


# Setup the plot 
f, (ax1, ax2) = plt.subplots(1, 2, figsize=(16,10)) 


# Rearrange and scale the data 
ax1.imshow(date1NDVI, cmap='brg',clim=(-@.1,0.7) ) 
ax2.imshow(date2NDVI, cmap='brg',clim=(-@.1,0.7) ) 


# Add some titles 
ax1.set_title('Autumn 2007 NDVI', fontsize=18) 
ax2.set_title(‘Autumn 2017 NDVI', fontsize=18) 


# Remove the axis labels for clarity 
ax1.yaxis.set_visible(False) 
ax2.yaxis.set_visible(False) 
ax1.xaxis.set_visible(False) 
ax2.xaxis.set_visible(False) 


Autumn 2007 NDVI 
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Now talking about this set of code for computing NDVI from the extracted 
subsets, np.seterr(divide=‘ignore’, invalid=‘ignore’) allows the devision by 
zero. Then we create a function makeNDVI to compute the NDVI from the 
red and NRI bands, ndvi = (b4.astype(float)-b3.astype(float)) / 
(b4.astype(float) + b3.astype(float)) is the formula to calculate the NDVI. 


ndvi[b4 == 32767] = np.nan function makes the NoData values to NaN so 
they get masked from the analysis, and then we return the function ndvi. 
date1NDVI & date2NDVI use the function to produce NDVI images, whereas 
f, (ax1, ax2) setups the plot. As explained before ax1.imshow & 
ax2.imshow rearranges and scales the data, whereas ax1.set_title & 
ax2.set_title adds title to our output image. 


set_visible(False) removes the axis labels for clarity, if we set it true then we 
can see the axis label as our previous output. 
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Autumn 2007 NDVI Autumn 2017 NDVI 


Now this is the clear output image, we can spot the NDVI change analysis in 
two images very easily. Now we have a look at how the values have changed 
in this 10 year period. To do this we simply subtract the two dates. Initially 
we plot a histogram to look at the spread of values. This is because 
differences and NDVI Are due to many things other than change. Primarily 
the largest difference here is due to climatic effects. It was much wetter at 
the start of this year than it was in the start of 2007 and so overall there is a 
general decrease in NDVI between these two time periods. However certain 
areas of change are obvious as well. There are more watering points 
(indicated by very low NDVI values) and there is more coal seam gas 
infrastructure. 


Now we can undertake a very simple normalisation by simply standardising 
each era so that all change is centred around zero. This makes the resulting 
change image slightly easier to interpret. 


ndviChange = date1NDVI - date2NDVI function computes the NDVI change 
between the dates, ndviChangeNorm function implements simple 
normalisation for the standardise change as we can se in figure on next 
page. 
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# Compute the NDVI change between the dates 
ndviChange = date1NDVI - date2NDVI 


# Standardise the change 
ndviChangeNorm = ((date1NDVI - np.nanmean(date1NDVI)) / np.nanstd(date1NDVI)) - \ 
((date2NDVI - np.nanmean(date2NDVI)) / np.nanstd(date2NDVI) ) 


# Setup the plots 
f, (axl, ax2) = plt.subplots(1, 2, figsize=(16,10)) 


# Rearrange and scale the data 
ax1.hist(ndviChange. ravel(), bins=256, range=(-@.5, @.5),color = "cyan") 
ax2.hist(ndviChangeNorm. ravel(), bins=256, range=(-3, 3),color = "blue") 


# Add some titles 
ax1.set_title('2017 - 2007 NDVI Change', fontsize=18) 
ax2.set_title('Standardised NDVI', fontsize=18) 


# Add some grids 
ax1.grid() 
ax2.grid() 


2017 - 2007 NDVI Change Standardised NDVI 
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As | explained before f,(ax1, ax2), setups the plot, ax1.hist & ax2.hist 
rearranges and scales the data, ax1.set_title & ax2.set_title adds title to our 
graph and ax1.grid() & ax2.grid() adds some grids to our plot. 


2017 - 2007 NDVI Change Standardised NDVI 
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Above graph clearly shows 2017 - 2007 NDVI change and other graphs tells 
us about standardised NDVI. 
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Now we can show the change in raster images, now we can make a image 
of the change in NDVI between the two dates. In a colour map we have 
selected below red indicates a decrease between the dates, green indicates 
an increase between the dates. 


plt.figure(figsize=(12,12) ) 
plt.imshow(ndviChangeNorm, cmap='RdY1Gn_r',clim=(-3.8,3.8)) 
# Add a title to the plot 


plt.title('1@ year Standardised NDVI Change', fontsize=18) 


Text(0.5, 1.0, '10 year Standardised NDVI Change' ) 


10 year Standardised NDVI Change 
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plt.figure(figsize=(12,12)) defines at size of figure we will receive as output, 
then in plit.imshow() we pass ndviChangeNorm which is change 
normalisation method as | explained before, cmap=‘RdYIGn_r’ which defines 
colour map and RdYIGn which means red, yellow and green is for the 
divergence maps and clim is used to set the limits of the current image. Then 
we will add a title to our plot using plt.title. 
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Above image is the final output for the change analysis performed, we can 
notice that aside from areas of infrastructure that have been put in and a 
small patch of clearing in the south-east corner there is an area of significant 
greening up in the north-east. 


4.2 Conclusion for Change Detection: 


Here we treated the data hosted on hittp://qid.auscover.org.au and use the 
RasterlO package to interact with the data and undertake some typical 
remote sensing tasks. In this part of our report we looked at changes 
between imagery collected at two different dates, there are many different 
ways of comparing the imagery typically you would want to run some sort of 
image classification approach which can be easily done but its beyond the 
scope of initial environment and workload. 


The Landsat data on TERN are pushed seasonally from DSITI and now cover 
all of Australia. They are mosaiced per state, and then a script on the server 
produces gdal virtual mosaics with national extent ( http:// 
qild.auscover.org.au/public/data/landsat/seasonal fractional cover/ 
ground cover/aus/ ). These are setup with internal web links so they work 
using GDAL's virtual driver system. What this means is that you can simply 
connect to the data on the server as if it were local, and only work with the 
bit you are interested in rather than downloading 10TB of data locally. 


For reading in a data subset from the remote file system, first we Compute 
the area we want to extract from the image and then we will pull the data 
subset from the server. In all these examples, subsets are extracted without 
having to read the entire file so you only pull the data you are working with 
back to your browser. In the case of this example, the (compressed) file is 
18GB in size but we just read the portion we want so the subset is available 
in under a second. 


Talking about Displaying a RGB image, Like any remote sensing package, 
we have to stretch the data into 8 bit space and set up an informative band 
combination.For more info on common band combinations (http:// 


landsat.usgs.gov/L8_band_ combos.php). 
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5. Ul Development 


UI ( User Interface ) development is the process of programming client-facing 
interfaces. The Ul development process includes writing code for images, 
animations, sliders, text fields, buttons, etc. 


Here for our Ul we use Kivy and KivyMd framework because this framework 
is easily compatible with python programming language. User interface is 
important to meet user expectations and support the effective functionality 
of your site. A well-executed user interface facilitates effective interaction 
between the user and the program, app or machine through contrasting 
visuals, clean design and responsiveness. The goal of a user interface is to 
maximize a user’s ability to access all necessary components of the intended 
application or website with less effort exerted by the user; as a developer 
you want to create something beautiful that somebody with any previous 
experience of interfaces would be capable of using. Using an interface 
should be an enjoyable experience, the user shouldn’t be stressed out or 
annoyed by UI, and its inability to present the necessary information in an 
aesthetically pleasing and efficient manner. 


When designing a UI for your site, it's important to consider the user's 
expectations in terms of accessibility, visual aesthetic and ease of use. An 
optimal mix of effective visuals and efficient responsiveness will improve 
your site's conversion rates, as it anticipates the needs of the user and then 
satisfies those needs. More specifically, here are the most important 
overarching elements of a great UI: 


1. Information Architecture 
2. Interactive Design 
3. Visual Design 


RESULTS 12 Reon oe sagt RESULTS 
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5.1 Ul Development Tutorial: 


from kivy.lang import Builder 

from kivy.properties import ObjectProperty 

from kivy.uix.screenmanager import ScreenManager, Screen 
from kivymd.app import MDApp 

from kivymd.uix.scrollview import MDScrollView 

from plyer import filechooser 

from training2 import CycloneSeverity 

from changedetection import ChangeDetection 
Builder.load_string(""" 

<NavigationDrawer> 


MDNavigationDrawerMenu: 


MDNavigationDrawerHeader: 
title: "Made In India" 
text; ©" 
source: "Tiranga.jpeg" 
spacing: "1@dp" 


MDNavigationDrawerDivider: 
MDNavigationDrawerItem: 


text: "Severity" 
icon: "“wind-power" 


on_press: 
root.nav_drawer.set_state("close") 
root.screen_manager.current = "scr 1" 


MDNavigationDrawerItem: 
text: "Change Detection" 
icon: "compare" 


on_press: 
root.nav_drawer.set_state("close") 
root.screen_manager.current = "scr 2" 


MDNavigationDrawerItem: 
text: "Settings" 
icon: "cog" 


on_press: 
root.nav_drawer.set_state("close") 
root.screen_manager.current = "scr 3" 


MDNavigationDrawerDivider: 


<WelcomeScreen>: 
MDScreen: 
MDBoxLayout: 
orientation: 'vertical' 
md_bg_color: (120/255.0,170/255.0,160/255.0,1) 


MDLabel: 
text: 'Welcome to UI Change Detection APP' 
halign: 'center' 
font_style:'H3' 
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To start Ul development using Kivy or KivyMd framework, first you have to 
install kivy in your working system using pip via your Terminal or Command 
Prompt. Just run “ pip install kivy ” or “ pip install kivymd ” on your 
terminal or cmd to install the following packages on your system and 
working directory. 


3 17 


In starting first we will import some important packages that we will need 
further. from kivy.lang import Builder which will import builder method, 
builder method is a helpful mechanism for using the builder pattern without 
writing boilerplate code and we can apply this annotation to a class or 
method. from kivy.properties import ObjectProperty, where 
ObjectProperty is a specialised sub class of the property class. from 
kivy.uix.screenmanager import ScreenManager, Screen imports screen 
manager and screens in our app. from kivymd.app import MDApp which 
inherits properties of a app from kivymd, from kivymd.uix.scrollview import 
MDScrollView we import scroll view which is configured to allow paging 
through views using swiping gestures by using the paging enabled props. 


Then from plyer import filechooser we import file chooser in our app 
through which we can choose a file to process from our PC. At last from 
training import CycloneSeverity & from changedetection import 
ChangeDetection we import logic from other py files, the two concepts | 
explained previously. Training.py is one of the py files from which | imported 
CycloneSeverity and changedetection is another py file from which | 
imported ChangeDetection. 


Now using Builder.load_string() we will start writing our basic code. First we 
will create Navigation Drawer using <NavigationDrawer> and under 
navigation drawer, we will get MDNavigationDrawerMenu which give a menu 
button in our app. Then under MdNavigationDrawerMenu we will put 
MDNavigationDrawerHeader which allow us make a header in navigation 
drawer and here under we will pass different attributes which md navigation 
drawer menu takes. As before under MDNavigationDrawerMenu we will add 
MDNavigationDrawerltem through which we can add multiple screen in our 
app and by attribute text we can add name to that screen, icon will add 
different icons you can get different icons from { https://fonts.google.com/ 
C O n Ss a 
selected=Material+Symbols+Outlined:menu:FILL@0:wght@400;GRAD@0:o0p 
$z@48&icon.platform=android&icon.query=application } this website. 


Attribute on_press defines what will happen when certain button gets 
clicked or here as we move towards different screen so it will provide us the 
functionality of different screens. 
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As here we give command root.screen_manager.current = ‘scr 1’ tells us 
that on pressing button severity in navigation centre opens screen 1 and 
screen manager manages these rotation of screens. Likely we will manage all 
3 screens using screen manager and provide on_press command to them so 
that they report as per their work. 


Then here we developed a welcome screen by using <WelcomeScreen> 
and define a MDScreen under welcome screen. Actually MDScreen is better 
version of screen so that’s the reason we use MDScreen instead of Screen. 
MDScreen will provide us here a whole new screen to work on. 


Under MDScreen we developed a MDBoxLayout which makes a box like 
structure in the screen. As our need we can provide specifications to our box 
layout like their background colour, size, position, etc. But here we just 
provided orientation: ‘vertical’ and md_bg_color: ( 120/255.0, 170/255.0, 
160/255.0, 1). Here you would be wondering what are those values written, 
in easy way here we can just provide name of some colour like ‘Red’, 
‘Yellow’, ‘Blue’, ‘Black’, ‘White’, ‘Green’, etc, you can study this in the 
documentation of KivyMD but here we used values where we assigned 4 
values in bracket. This setup represents RGB colours first value is for Red, 
second value is for Green, third value is for Blue, and fourth value is for 
Opacity which defines opacity of background colour. 


Benefit of using this setup is we can mix these colours to make any colour of 
our choice, if here we put the value (1,0,0,1) then background will be Red, if 
we put the value (0,1,0,1) then output will be green and background colour 
will be Green. The value we provided will provide a greenish colour with the 
pinch of blue colour in it which | will show you later on. 


Then under MDBoxLayout we develop a MDLabel, you must be thinking 
why MDLabel is under MDBoxLayout, its because we want to write 
something under our Box Layout so we use labels here. 


Here we provide 3 specifications inside MDLabel, first is text: ‘Welcome to 
UI Change Detection App’ which will come writing on screen as soon as we 
open the app. Then second is halign: ‘center’ which aligns text to the center 
of the screen. Third is we provide font_style: ‘H3’, which converts text to 
font heading type 3 from default. 
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MDFloatingActionButton: 
icon: 'arrow-right' 
pos_hint:{'center_x': 0.9, 'center_y':@.2} 


on_press: 
root.manager.transition.direction = 'left' 
root.manager.current = 'settings' 
<MainScreen> 
MDScreen: 
MDTopAppBar: 
id: top 


title: app.title 
pos hints {"top": i) 


elevation: 


4 


left_action_items: [["menu", lambda x: 
nav_drawer.set_state("open")]] 


MDNavigationLayout: 


MDScreenManager: 
id: screen_manager 


MDScreen: 
name: "scr 1" 


MDLabel: 


text: "Cyclones are caused by atmospheric 
disturbances around a low-pressure area 
distinguished by swift and often destructive air 
circulation. Cyclones are usually accompanied by 
violent storms and bad weather. The air circulates 
inward in an anticlockwise direction in the 
Northern hemisphere and clockwise in the Southern 
hemisphere. Cyclones are classified into five 
different levels on the basis of wind speed. They 
are further divided into the following categories 
according to their capacity to cause damage:-" 

halign: "center" 

pos_hint: {"center_x":0.5, "center_y":0.85} 

bold: True 

font_style: "H6" 


Image: 


source: "Untitled.jpeg" 
pos_hint: {"center_x":@.5, "center_y":@.7} 


MDFillRoundFlatButton: 


text: 'Choose File' 

pos: root.width*@.15, root.height*@.44 
md_bg_color: "black" 

text_color: 'white' 

on_release: root.file_chooser() 


MDLabel: 
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Then after MDLabel we develop a button. There are varieties of button which 
you can find reaching to this link { https://kivymd.readthedocs.io/en/1.1.1/ 
components/button/ }, here we will be using MDFloatingActionButton to 
change the welcome page. In this button we can add icons too with same 
process as previous, so here we added 3 arguments in button first is icon: 
‘arrow-right’, second is pos_hint: {‘center_x’:0.9, ‘center_y’:0.2} through 
which we buttons position on screen and third is on_press: 
root.manager.current = ‘settings’. 


Now from this point we start to develop our Main Screen by passing 
argument <MainScreen> and adding widgets to our main screen. Under 
main screen first we define a screen using MDScreen and under screen we 
define MDTopAppBar which develops a bar at the top of the screen. Under 
top app bar we pass 5 arguments first is id: top, second is title: app.title, 
third is pos_hint: {“top”:1}, third is elevation: 4 and last one is 
left_action_items: [["menu", lambda x: nav_drawer.set_state(“open")]] 
which defines what action item menu button will contain in the top app bar 
where we define nav_drawer.set_state(“open”). 


Then after MDTopAppBar we define MDNavigationLayout under navigation 
layout we define MDScreenManager with a single argument id: 
screen_manager. Then under screen manager we define screen as 
MDScreen and name it as scr 1 and then under screen we define MDLabel. 
We provide a text argument under label and pass the text that we wanted to 
wright, after we give pos_hint which specifies location of text and also 
halign argument which keep text in middle of the screen. At last we provide 
font_style as H6. 


After MDLabel we define Image through which we add a image to our 
screen. Under image we pass 2 argument where first is source in which we 
provide name of the file and then we provide pos_hint which defines position 
of the image in the screen. 


Then we develop a button using MDFillRoundFlatButton through which we 
will choose files from our PC. Under button we pass 5 arguments first is text 
that appears as the name of button, second is pos which defines where 
button is located on the screen, third is md_bg_color which defines colour of 
button , fourth is text_color which defines the colour of text on the button 
and at last we define on_release: root.file_chooser() which defines what 
button will execute which we will discuss afterwards. 
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id: selected_path 

texes oe 

pos_hint: {"center_x":@.57, "center_y":0.5} 
theme_text_color: "Custom" 

textcolor:. (1) 0; 8; L) 

bold: True 


MDRectangleFlatIconButton: 
icon: 'wind-power' 
text: 'Check Severity of Cyclone' 
pos: root.width*@.15, root.height*@.37 
line_color: '‘'red' 
on_release: root.cyclone_severity() 


MDScreen: 
name: "scr 2" 


MDLabel: 
text: "Change detection for GIS (geographical 
information systems) is a process that measures 
how the attributes of a particular area have 
changed between two or more time periods. Change 
detection often involves comparing aerial 
photographs or satellite imagery of the area taken 
at different times. The process is most frequently 
associated with environmental monitoring, natural 
resource management, or measuring urban 
development. Change detection is done by comparing 
two raster images that were taken at different 
times but which cover the same area. As the images 
cover the same area, the images overlay each 
other. Imagine two grids stacked on top of each 
other. It is then a matter of comparing whether 
the value of a pixel in the new raster is the same 
as the value of the pixel in the old raster. 
Pixels that have changed are then marked. The 
output is usually a raster that covers the same 
extents as the two images with the changed areas 
highlighted." 
halign: "center" 
pos_hint: {"center_x":@.5, "center_y":@.85} 
bold: True 
font_style: "Hé6" 


Image: 
source: "uoo5z.png" 
pos_hint: {"center_x":@.2, "center_y":@.62} 


MDFillRoundFlatButton: 
text: 'Choose File 1' 
pos: root.width*@.15, root.height*@.34 
md_bg_color: "black" 
text_color: 'white' 
on_release: root.file_chooser() 
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MDLabel: 


id: selected_path 
texte: 23 


pos_hint: {"center_x":@.57, "center_y" 


theme_text_color: "Custom" 
text color: 2} 6). Oy 1) 
bold: True 


MDFillRoundFlatButton: 


text: 'Choose File 2' 


pos: root.width*@.15, root.height*@.18 


md_bg_color: "black" 
text_color: ‘white' 
on_release: root.file_chooser() 


MDLabel: 


id: selected_path 
text: ** 


pos_hint: {"center_x":0@.57, "center_y" 


theme_text_color: "Custom" 
text color: (1; 0; @;.-2) 
bold: True 


MDRectangleFlatIconButton: 


icon: 'compare' 
text: 'Detect Change' 


pos: root.width*@.35, root.height*®.26 


line_color: 'red' 
on_release: root.change_detection() 


MDScreen: 
name: "scr 3" 


MDLabel: 


text: "Settings" 
halign: "center" 
font_style: "H4" 


MDNavigationDrawer: 
id: nav_drawer 


radius: 


(85; 16;. 16, 6) 


NavigationDrawer: 
screen_manager: screen_manager 
nav_drawer: nav_drawer 


class WelcomeScreen(Screen): 


pass 


class MainScreen(Screen): 


def file_chooser(self): 


filechooser.open_file(on_selection = self.selected) 
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20.5} 


20.5} 


After button we define MDLabel again in which we pass 6 arguments where 
first is id: selected_path, second is text, third is pos_hint, fourth is text_color, 
fifth is theme_text_color and last is bold. When we press the choose file 
button on the screen and select certain file from our pc then this label will 
display the path of the file. 


And then under screen 1 we process our last widget which is 
MDRectangleFlaticonButton — which will process our cyclone severity 
algorithm which we learnt previously in our report. Under button we define its 
specification and at last we define on_release: root.cyclone_severity which 
we will understand later in this report. 


Now we will define our second screen in main screen using MDScreen with 
name scr 2. Then as we did earlier we will add label to our and provide text 
we want to display on screen, provide there positions hint, font_style, etc 
arguments under MDLabel. 


Then we will again add Image as we did before and give source of image 
and its position hint. After that we will define a button 
MDFillRoundFlatButton to choose file from pc and provide its different 
needed specifications under it. After that we define MDLabel which will print 
the path of file on screen which we will choose from pc. 


And as same again we will add a MDFillRoundFlatButton to choose file 
because here we are going to show change detection in raster images So we 
will need 2 different button to choose 2 different files and as same again we 
add MDLabel after button which will print the path of file. 


Now we will add one more button to scr 2 which will process both the 
images and show change detection both the images. We chooses 
MDRectangleFlatlconButton for this process and specify on_release: 
root.change_detection(). That’s how this button will process our change 
detection file. And development of screen 2 is completed. 


Now we will develop our last screen which is screen 3 which will dispaly app 
settings window. So again we add MDScreen and name it as scr 3. After this 
we define MDNavigationDrawer and gives it id as nav_drawer and defined 
its radius. 


Now its time to end our Builder method and pass every logic that we 
mentioned so first we passed WelcomeScreen. Then made a class 
MainScreen and defined file_chooser as self inside the class. 
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def 


def 


def 


pass 


selected(self, selection): 
if selection: 
self.selected_path = self.selected 
self.root.ids.selected_path.text = selection[@] 
pass 
cyclone_severity(self): 
try: 
if hasattr(self, 'selected_file'): 
file_content = CycloneSeverity(self.selected_path) 
print("File content:") 
print(file_content) 
else: 
print("No file selected.") 
except Exception as e: 
print(f"Error processing the file: {e}") 


CycloneSeverity.open_file(on_selection = self.selected_path) 


change_detection(self): 
try: 
if hasattr(self, 'selected_file'): 
file_content = ChangeDetection(self.selected_path) 
print("File content:") 
print (file_content) 
else: 
print("No file selected.") 
except Exception as e: 
print(f"Error processing the file: {e}") 


ChangeDetection.open_file(on_selection = self.selected_path) 


class NavigationDrawer(MDScrollView) : 
screen_manager = ObjectProperty() 


nav_ 


drawer = ObjectProperty() 


class App(MDApp): 


def 


def 


__init__(self, **kwargs): 
self.title = "Change Detection UI" 
super().__init__(**kwargs) 


build(self): 
self.theme_cls.primary_palette = "Green" 
self.theme_cls.theme_style = "Light" 

sm = ScreenManager() 

sm.add_widget (WelcomeScreen(name='menu' ) ) 
sm.add_widget (MainScreen(name='settings' ) ) 
return sm 


if _name__ == "__main__": 
App().run() 
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Then inside file_chooser we pass an argument 
filechooser.open_file(on_selection = self.selected) which defines that as 
we click the file chooser button file section screen will appear on its own. 


And in next step we defined selected with arguments self and selection and 
then inside selected we define a if loop which says that self.selected_path = 
self.selected which means selected path is selected. And then this line 
self.root.ids.selected_path.text = selection[0] defines that the path of the 
selected path will be displayed as text on screen which | explained you while 
explaining builder method. 


Then we define cyclone_severity as self and under that we define try with if 
else loop. You can see that we defined is with has attributed and based self 
and selected_file in it and then we gave logic of our code will process the 
given path in certain file which is file_content = 
CycloneSeverity(self.selected_path) and then we print the content of the 
file. And then in else we passed an exception e which will print the error 
message if no files are selected. 


As same as cyclone_severity we defined change_detection as self and as 
same as above we defined if else loop inside try and defined an exception as 
e to print the error message if file not selected and then passed the class 
MainScreen. 


After that we again made a class NavigationDrawer and passed 
MDScrollview in it so that it can inherit the property of scroll view. Under 
navigation drawer we passed two functions screen_manager and 
nav_drawer and defined it as ObjectProperty(). 


At made a class App and passed it with MDApp, this main class which is 
important to execute the app. Under it we defined __init__ in which we gave 
our app a name. 


And then we defined our build method with self were we theme class primary 
palette as green and theme style as light. Then we defined our screen 
manager as sm, and we added our both the screens WelcomeScreen and 
MainScreen to screen manager with method sm.add_widget. then we return 
the sm and call method App().run() which runs our app successfully. 
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@ python = Window ® Sat 5 Aug 5:41PM 1 


Change Detection UI 


Welcome to Ul Change Detection APP 


This is our running app, you can see welcome page appears as soon as 
we open the app. 


Cyclones are caused by atmospheric disturbances around a low-pressure area distinguished by swift and often destructive air circulation. Cyclones are usually accompanied by violent storms 
{and bad weather. The air circulates inward in an anticlockwise direction in the Northern hemisphere and clockwise in the Southern hemisphere. Cyclones are classified into five different levels on 
the basis of wind speed. They are further divided into the following categories according to their capacity to cause damage:- ) 


Wind Speed in Km/h Damage Capacity 


linima 


Choose File 


=A Check Severity of Cyclone 


This is our first screen of severity, where we can find severity of 
cyclone. As soon as we choose the file and press button ‘check severity 
of cyclone’ then result will process on right side of screen where we can 
see the image of cyclone and stage of cyclone would be written with 
prediction accuracy. 


Made In India 


={ Severity 


A) Change Detection 


2% Settings Wind Speed in Km/h 


Here we can see the menu of our app which has 3 different fields, first 
Severity, second Change Detection and third Settings. 
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Change detection for GIS (geographical information systems) is a process that measures how the attributes of a particular area have changed between two or more time periods. Change 
detection often involves comparing aerial photographs or satellite imagery of the area taken at different times. The process is most frequently associated with environmental monitoring, natural 
resource management, or measuring urban development. Change detection is done by comparing two raster images that were taken at different times but which cover the same area. As the 
images cover the same area, the images overlay each other. Imagine two grids stacked on top of each other. It is then a matter of comparing whether the value of a pixel in the new raster is the 
| same as the value of the pixel in the old raster. Pixels that have changed are then marked. The output is usually a raster that covers the same extents as the two images with the changed areas. 
highlighted. 


Image Differencing 


Choose File 1 


A) Detect Change 


Choose File 2 


This is our second screen of Change Detection, where we can find 
change detection in raster images. As soon as we choose the file and 
press button ‘Detect Change’ then result will process on right side of 
screen where we can see the image which shows the change image 
with red, yellow and blue shades on the image. 


So here we end UI Development Tutorial. 
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6. Advantages of using CNN for Machine 
Learning and Deep Learning 


At the end | would like to explain you some major benefits of CNN in field of 
machine learning and deep learning. Deep learning is a form of machine 
learning that requires a neural network with a minimum of three layers. 
Networks with multiple layers are more accurate than single-layer networks. 
Deep learning applications often use CNNs or RNNs (recurrent neural 
networks). 


The CNN architecture is especially useful for image recognition and image 
classification, as well as other computer vision tasks because they can 
process large amounts of data and produce highly accurate predictions. 
CNNs can learn the features of an object through multiple iterations, 
eliminating the need for manual feature engineering tasks like feature 
extraction. 


It is possible to retrain a CNN for a new recognition task or build a new 
model based on an existing network with trained weights. This is Known as 
transfer learning. This enables ML model developers to apply CNNs to 
different use cases without starting from scratch. 

1. Human supervision not required 

2. Automatic feature extraction 

3. Highly accurate at image recognition & classification 

4. Weight sharing 

5. Minimizes computation 

6. Uses same knowledge across all image locations. 


7. Ability to handle large datasets 


8. Hierarchical learning 


These are some of the most significant advantages of convolutional 
neural networks (CNNs). 
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7. Conclusion 


In this project and report, we came across a variety of concepts that are part 
of the Convolutional Neural Network. Each part has its own importance and 
combined they make up a highly efficient and good performance neural 
network known as CNN. They are an essential part of deep learning 
applications. 


CNN is computationally efficient, It performs parameter sharing and uses 
special convolution and pooling algorithms. CNN models may now run on 
any device, making them globally appealing, It finds the relevant features 
without the need for human intervention. It can be utilized in a variety of 
industries to execute key tasks such as facial recognition, document 
analysis, climate comprehension, image recognition, and item identification, 
among others. By feeding your data on each level and tuning the CNN a little 
for a specific purpose, you can extract valuable features from an already 
trained CNN with its taught weights. 


Convolutional Neural Net is a popular deep learning technique for current 
visual recognition tasks. Like all deep learning techniques, CNN is very 
dependent on the size and quality of the training data. Given a well prepared 
dataset, CNNs are capable of surpassing humans at visual recognition tasks. 
However, they are still not robust to visual artifacts such as glare and noise, 
which humans are able to cope. The theory of CNN is still being developed 
and researchers are working to endow it with properties such as active 
attention and online memory, allowing CNNs to evaluate new items that are 
vastly different from what they were trained on. This better emulates the 
mammalian visual system, thus moving towards a smarter artificial visual 
recognition system. 
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